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, The structure of the isovector dipole resonance in neutron-rich calcium isotope, 

""Ca, has been investigated by implementing a careful treatment of the differences 
• of neutron and proton radii in the continuum random phase approximation (RPA). 

■ The calculations have taken into account the current estimates of the neutron skin. 
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The estimates of the escape widths for direct neutron decay from the pygmy dipole 
resonance (PDR) were shown rather wide, implicating a strong coupling to the 



OO 

■ continuum. The width of the giant dipole resonance [GDR] was evaluated, bringing 

00 , 

<^ ', on a detailed discussion about its microscopic structure. 
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PACS numbers: 21.10.Pc, 21.60.-n, 24.30.Cz, 24.30.Gd 
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^ ■ I. INTRODUCTION 

> 

^ \ The investigation of the nuclei lying far from the j3~ stability line has been an interesting 

c3 ■ and active field of the nuclear physics in the last two decades. In that region, a great 
number of exotic features are observed like halo/skin formation, intruders levels, new magic 
numbers and new kinds of collective excitation, the so-called soft and pygmy resonances. 



These observations have 
the j3~ stability valley 



breed the review of successful theoretical tools in nuclei around 



mm 



The description of the microscopic structure of the exotic nuclei is a current topic of study 
and several recent works had the concern of describing the giant resonance {GR) in neutron- 
rich nuclei. The main questions treated are the giant dipole resonance {GDR) behavior 
and the appearance of the pygmy dipole resonance {PDR) in a nucleus with large ratio of 
neutron to proton number {N/Z). The PDR appears in medium and heavy neutron-rich 
nuclei, and, within the hydrodynamic sense, they are probably due to oscillation of the 
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neutron excess against the core (A^ = Z nucleus). The effects of the neutron excess in dipole 
resonances, in neutron-rich nuclei, have been extensively studied in the literature by many 

group. QBBBayyQQ , but a complete understanding about the resonance 

structure has not been reached up to now. Some questions about the nature of GDR and 
PDR are still kept open. The preservation of the structure of the GDR as the neutron 
number increases, the com pos ition of its decay modes and the need to go beyond Ip — Ih 



configuration to explain it 



ll| are some of these subjects of theoretical interest. In addition 



to these interesting problems on nuclear structure, the PDR in neutron-rich nuclei plays 
an important subject in the r-process nucleosynthesis [l^ Q. In this topic, the dipole 
strength and the widths of the PDR and GDR are important ingredients to understand the 
radiative neutron capture, because they are directly related with the competition between 
the direct and statistical mechanisms in the neutron capture process. In order to perform 
the studies on neutron-rich nuclei, we choose to analyze the isovector dipole resonance in 
the neutron- rich calcium isotope ^^Ca. This nucleus was not reached experimentally yet, 
but its microscopic structure has been investigated recently in different versions of RPA 
calculations flE.SS, relativistic RPA (RRPA) Q and phonon damping model (PDM) 

In this work we are interested in studying the role of continuum, as well as the decay 
width, in the pygmy and giant dipole resonances in ^^Ca nucleus. In these extremely fragile 
nuclear systems, close to the neutron drip line, the neutron density extends far away from the 
proton density, forming the 'neutron skin'. Thus, to perform these continuum calculations, 
we have implemented a careful treatment of the differences of neutron and proton radii in a 
version of the continuum RPA approach of the Ref. Q|. 

In the Sect. II of this paper, we describe the theoretical approach used in the calculations. 
In the Sect. Ill we present and comment our results. 



II. THEORETICAL APPROACH 



The continuum effects in our microscopic calculations are taken into account through a 
discrete particle-hole basis which accomodates the single-particle resonance widths, resulting 
in a diagonalization of RPA-like complex matrixes of standard size Q|. The excited states 
are implemented in the particle-hole excitation space, and the open channels correspond to 
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unbound particle-hole states with complex energies, of which the imaginary parts give the 
single-particle escape widths. We assume the nuclear hamiltonian of the form H = Hq + 
\ ^ijtj ^ij with a mean field part Hq and a residual two-body force so that {E^, — H) | z/) = 
is satisfied for the nuclear excited states | v) with energy E,^. Using the orthogonal and 
complementary projectors Q and P, the particle-hole space can be splited onto normalized 
bound (or unbound) particle-hole states and residual continuum states, respectively. Then, 
using the projection operator formalism, we have the equations: 

[E^ - Hqq - Hqp I Hpq]Q I u) = Hqp I x^) , (1) 

where \x^) is the appropriate scattering solution of Hpp: 

[Eu - Hpp] I x^) = 0. (2) 

The relevant continuum effects in P space is taken by accounting the single-particle resonance 
in unbound particle-hole states (taken out for the Q space) that is reached making the 
approach of ignoring the two body force in P space. In this way, the Eq.Q is a single- 
particle equation and the continuum self-eriergy {Hqqp ^ -i-i^ijjppp -^opq) dresses the single- 
particle resonance with escape effects Thus, we define the complex particle-hole 
modes | -R„) with complex energies = 5„ — | i Fj^ which account for continuum escape 
effects: 

- Hqq] \Rn) = and [e: - H^q] \ P„) = , (3) 

Q^^^^j^ \Rn){Rn\HQp\x^) 

where Hqq = Hqq + Hqp ^ +i'^-Hpp ^PQ ^^"^ complex modes satisfy the orthonorma- 
lization relation: {R^i \ Rn) = The matrix element in Eq.Q is an escape amplitude 

T \ I R \f \ o\ I ^ 

related to the imaginary part of The strength function Sf{E) = Y,n , "A ^ ' t , ^2 
calculated by making the approach that the excited states | u) can be well described by the 
Q I u) component. In the particle-hole matrix element calculation of 1-body operator Fx of 
the strength function Sf{E), we have assumed the form 

FjM = e^YjM (5) 
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where is the nucleon effective charge. For isovector dipole transition, we have ej,(^) = — ^ 



(f)- 

The complex particle-hole modes | in Eq.(j2I) are solved by a diagonalization of the 
discrete RPA equations in Q space: 

' Mf""l^J-] (6) 

y-B -A J yv'' J V / 

where 

Aphp'h' = i^p — £h)Spp'5hh' + Vph'hp' ; Bphp'h' = Vpp'hh' (7) 

and 6p {sp = Bp — \ i Tp) are the complex energies of the single-particle resonances 
Q]. Diagonalizing the complex Eq.®, we have the complex eigenvectors | -R„) (given 
by complex X^^ and 1^*^ amplitudes) and complex eigenvalues e„. Thus, in this discrete 
particle-hole subspace, the escape width is associated with the contribution of all allowed 
unbound particles coupled to their respective single-holes. The partial escape width for each 
single-hole can be approximated for ~ Y.p X^h ■> what gives a good estimate of the 
escape width composition for the population of several single-holes in the residual nucleus. 

The discrete single-particle energies are evaluated by solving the Schrodinger equation 
with Woods-Saxon potential, including the centrifugal and Coulomb (as a uniformly charged 
sphere) terms. The positive single-particle energy and its respective width are calculated 
in a projection technique to continuum discretization approach jisl . The potential 

parameters of ^''Ca were adjusted following the systematics of nucleon binding energy for 
this nucleus, because there is no experimental data available for this nucleus. The '^^Ca 
is expected to have a small, but positive, neutron energy separation and a large proton 
energy separation {Sn ~ 3.5 MeV and Sp ~ 25 MeV). As a consequence, the single-particle 
potential of proton should be deeper than the neutron one. Moreover, Since 60Ca is already 
close to the neutron drip line, the next nucleus with full neutron subshell, ^°Ca, should be 
not stable against neutron emission. In TABLE the potential parameters are displayed. 
The single-particle energies calculated with these parameters are very similar with those 
resulting of self-consistent Hartree-Fock (HF) calculation with SHI Skyrme interaction 

nn 

|12l,llSl|, except for the first neutron excited state, namely lgg/2 level, which was evaluated 
to be unbound {sigg/^ — 0-9 MeV), although, in other calculations, it has been foreseen in a 
negative energy close to the border of the potential well. However, these differences do not 
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TABLE I: The Woods-Saxon central and spin-orbit parameters used in the calculations. 



VR{r) = Vmf{r) ; Vis{r) = V^i, 




\r{r)l ■ s 




; Vb/s = 


6.54 MeV 


Particle VoFiiMeV) 


i?(/m) 




V 54.5 


4.30 




TT 65.5 


4.62 





affect our results for the calculated escape widths, because this is a very narrow resonance 
(Fi,,^, < 1.0 keV). 

The RPA calculation is done by utilizing the Landau-Migdal residual interaction: 

V,h{n, n) = Coif in) + f\n)Tl ■ rf + 

+^-^igin)+ g'{n)T^ ■T^)]5{¥t (8) 
where /, /', g and g' are dimensionless and density dependent parameters: 

F(r) = F"^ + (F*" - F'^^)e(r) . (9) 



The set of the interaction parameters is similar to that used in Ref. jl^l, |20|, and it was 
adjusted to eliminate the spurious state 1~ (see TABLE lU]) and to reproduce the first 3~ 
excited state around 2.0 MeV. 

We have considered a two-parameter Fermi distribution to represent the density depen- 
dence ^(r) in Eq.(jni). In previous calculations, we had assumed: 

= 1 + Jr-nya ' (10) 

where R and a are half-density radius and diffuseness, respectively. Nevertheless, in the 
present context, it is more appropriate to separate this function into neutron and proton 
parts: 

= ^Ur) + jUr) , (11) 

where each part is given by: 

Here k = z/(7r) for neutron (proton). The Fermi distribution parameters are adjusted to 
reproduce the neutron and proton expected root mean square radii for the nucleus consid- 



ered. It has been observed that the neutron root mean square radius {J{r'^)v) becomes 
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TABLE II: Parameters of harmonic oscillator radial wavefunctions and residual interaction. 







60Ca 






v/(r^).(/m) R,{fm) 


auifm 


) V{r%ifm) 


R-n{fm) aT,{fm) 




4.03 4.34 


0.60 


3.62 


3.68 0.60 


3.90 


K{fm) 






Kifm) 




2.09 






2.11 




Co{MeV fm^) 


fex 




g 


9' 


300.0 -0.002 


-2.1 


0.76 


2.3 0.51 


0.70 



larger than the respective proton radius {yir'^)^) as the neutron number increases, keeping 
constant the proton number. The difference between those radii, Ar^,, = ^ {r'^)u — \J {r'^)TT, 
gives the measurement of the 'neutron skin'. Some estimates of Ar^^^ have been extracted 



from experimental charged radii 



21 



22 



24| . We perform the Rk and adjustment 



using the fact that for a two-parameter Fermi distribution {pk{f^) oc ^fc(r)) the mean square 
radius is given by j^]: 

j^Pk{r)r^dr 3 . 7 



(13) 



!^ Pk{r)r^dr 5"~" ' 5 

In our RPA calculation, we represent the radial single-particle orbits by harmonic oscil- 
lator radial wave functions Ri{r) which are characterized by a size parameter 6, and we have 
adopted a two-parameters Fermi density dependence in the residual interaction. Since the 
neutron and proton densities are too different in exotic nuclei, we must pay special atten- 
tion with the parameter h. Therefore, we have performed a careful choice of the parameter 
h to analyze the microscopic structure of this exotic nucleus. In this way, using the radial 
distribution of harmonic oscillator, pfc(r) = PkJ2ei<ep (2ji + 1) -Rf ^fc), in order to calcu- 
late (r^)fe, we satisfy this requisition by assuming different size parameters for neutrons and 
protons. 



bl 



(3) 



where Xk = N {Z) for k = and (r^)^ is given by Eq.dlSI). In the TABLEjlT 
displayed the Rk and adjustment based on some recent radii estimates 



(14) 



we have 
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FIG. 1: The calculated strength function Sf{E) by continuum RPA approach for the El excitation 
in ^^Ca nucleus. The wide peak in ~ 8.6 MeV is in the energy interval where the PDR is expected. 

III. DISCUSSION AND CONCLUSIONS 



In this section, we present and discuss some results which we have obtained by using 
lp~lh continuum RPA approach, as described in the previous section, for isovector dipole 
electric transition for ^^Ca. The single-particle (hole) energy levels were taken from the lsi/2 
hole up to lg7/2 particle orbits. In this case, the core (A^ = Z = 20) is formed by neutrons 
and protons filling the energy levels up to sd shell and the neutrons excess occupying the pf 
shell above the core. The calculated strength function Sf{E) is displayed in the FIG.HJ 

In the ^°Ca nucleus, the GDR is located around an energy (~ 20 MeV ) above the proton 
{Sp ~ 8 MeV) and neutron separation energy {Sn ~ 16 MeV) HQ- 

Consequently, the 

proton and neutron emission competes by decay of this mode. Taking into account the 
neutron-rich nucleus (®°Ca), the theoretical predictions have displayed that PDR appears 
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TABLE III: The evaluated escape widths, and the estimates of the partial escape widths for each 
single-hole, of the three main peaks that compose the PDR in ~ 8.6 MeV of excitation energy 
(see FIG. HJ exhausting about 8% of the EWSR. 

En -\iTl{MeV) single - hole (MeF) 
8.05 -i 1.29 (1/5/2)1^ 2.26 

(2pi/2)^ 0.39 

(2p3/2). 0.02 

8.69 -i 1.35 (1/5/2)1. 0.05 

(2^1/2),. 2.65 

(2p3/2). 0-20 

8.94 -i 0.40 (1/5/2)^ 0.04 

(2^1/2),. 0.09 
(2p3/2). 0.71 



below the GDR as a result of neutron excess. The PDR and GDR are expected to appear 
below 10 MeV and 15 - 20 MeV 3,0,330,1111, respectively. The GDR is shifted down 
due to the increase of the mass. However, its microscopic structure is found to be composed 
mainly by the single-holes that belong to the most internal structure of the core, showing 
that they belong to the same category of excitations of the ^^Ga. As we have considered 
Sn ~ 3.5 MeV and Sp ^ 25 MeV for the ^^Ga nucleus, the one-neutron channel is presumed 
to be open in the PDR region. 

Acccd,,, .0 p.ev,ous RPA calculate. 0038, ou. ca.c„.a«o„. also p.ed.ct a con- 
siderable strength in the energy region below 10 MeV , what can be observed by the presence 
of the broad neutron peak in low-lying energy in FIG. H The broad width is due to the 
fact that the energy of the resonance to be above of the small neutron separation energy, 
implicating a strong coupling of the external neutrons to the continuum region. These 
low-lying energy states are the natural candidates for PDR because they have a dominant 
contribution of 'neutron skin' {pf shell), in agreement with hydrodynamic interpretation. 
In FIG. Qthe wide peak in ~ 8.6 MeV is composed by the overlap of three main peaks 
(see TABLE Uni) exhausting about 8% of the Energy Weighted Sum Rule (EWSR), which 
are composed by transitions involving neutron of 'skin' {2p and If 3s, 2d and Ig). This 
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TABLE IV: Estimates for the mean values of the escape widths for the PDR. 
En -\iv\{MeV) single - hole T^^MeV) 
8.6 -i 1.3 (1/5/2). 0.7 

(2^1/2);. 1.6 
(2p3/2). 0.2 



aspect becomes more clear when the structure of the RPA wave function of these peaks is 
observed. The wave functions of the states at 8.05 MeV , 8.69 MeV and 8.94 MeV are com- 
posed by the main neutron Ip — l/i configurations: (76% 1/5^22*^3/2) +14% 2p^^^3si/2j ), 
(93% 12^-/22^3/2)^ + 6% 12^-/22^3/2) J and (88% |2p-/224/2)^ + 9% 1 1/7-/^21^9/2) J, respec- 
tively, showing that each one of these peaks has a dominance of two Ip — Ih configurations. 
These results are in agreement with the HF + RRPA calculation of the Ref . |lo| , where 
they have noticed that the one, or at most two neutron Ip — Ih configurations, determine 
the structure of the low-energy states for isovector dipole resonance in ®°Ca, in contrast to 
the structure of the collective states which are characterized by a coherent superposition of 
many Ip — Ih excitations. However, the PDR is composed by the overlap of those three 
main states, and the collectivity can be understood through the participation of all particle- 
hole components that belong to this overlap. The estimates of the partial widths of 
the main neutrons single-holes that populate the PDR are presented in TABLE IIIII We 
have also evaluated an average width for the PDR, calculated as the weighted average of 
the widths of the peaks that compose it (see TABLE llVjl . the weights being the intensities 
of each peak relative to the 8% of the EWSR that they exhaust. In comparision to these 
predictions to the PDR, we can note an apparent lack of neutron escape in the region of 
energy around the 15 MeV , namely GDR region, where the peaks are very narrow (see 
FIG.^). These peaks are constituted principally by bound Ip — Ih configurations of protons 
and their narrow widths are due to the small neutrons amplitudes in this energy interval, 
refiecting a unexpected suppression of the neutron emission. This result is mainly due to 
the change of status of the neutrons single-particle levels, because those neutrons of pf shell, 
that belonged to the configurations of the particle states in ^°Ca nucleus, have become hole 
states in the neutron-rich nucleus, and their contributions for the GDR of the ^'^Ca are small 
(they are mainly responsible for the PDR appearance). Furthermore, the composition of 
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Energy (MeV) 



FIG. 2: Comparison among our results (solid line) with the one of other works: HF + RPA 
calculations by Catara et. al. |^ (dotted line), Hamamoto et. al ^ (dashed line) and HF + RRPA 
calculations by Vretenar et. al. (crosses). Our calculation is performed by including an 

arbitrary constant width T = 1.0 MeV, instead of calculated single-particle widths used in the 
calculation presented in FIG. ^ 

these narrow peaks is owed mainly to the bound Ip — lh pairs of protons. This neutron lack 
can indicate that the GDR decay should also have more complex contributions than Ip— Ih 
excitations. The analysis of the underlying structure of the resonance in the exotic nuclei 
may help us understand such disagreements because the GDR could have a more complex 
structure than Ip — Ih configurations, what could contribute to the spreading width. To 
simulate these more complex structures of the GDR, we display in the FIG.|21our calcula- 
tion using a constant width F = 1.0 MeV, instead of the single-particle widths used in the 

calculation presented in the FIG. [Tl and we compare the results with the one of the Refs. 
1 

101 . Moreover, the two-neutron separation energy (S'2n ~ 7 MeV) is expected to be too 
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smaller than Sp in the neutron rich nuclei 



27|. Then the two- neutron channel is presumed 



to be open above the PDR and below the GDR energy, and it is expected to constitute 
an important part of GDR decay. As our continuum RPA calculations have taken into 
account only Ip — lh configurations, this two-neutron emission is not considered. Therefore, 
the inclusion of 2p — 2h configuration in the calculations would enhance the neutron emission 
by the GDR. 

At present, it may be difficult, or else impossible, to make the experimental measures of 
these widths, but it could be worthwhile doing a theoretical effort for a better understanding 
of the microscopic structure of these exotic nuclei, that would allow us to accomplish some 
estimates. In Ref. [3l the low energy dipole strength, in capture cross section calculation, 
is folded by a Lorentzian curve of width Tpdr = ^gdr {Epdr/ Egdr)'^ , and it is added to 
the GDR strength to satisfy the classical sum rule for the total El strength. We could use 
this relation to perform an evaluation of the Tgdr- In TABLE llVI we show the mean values 
of the energy (8.6 MeV) and escape width (2.6 MeV) for the PDR. The large value of the 
escape width could suggest that it is the main component of its total width, or Tpdr = 2.6 
MeV. On the other hand, our calculation gives a narrow GDR centered in 15 MeV. In 
spite of our calculation does not contain 2p-2h or more complex excitations, we can use these 
results for PDR and GDR, together the above expression for the TpDp, to evaluate 
the width Tgdr- Proceeding this way, we obtain Tgdr = 7.9MeV. This result shows that 
the GDR is a very wide resonance, but with a small Ip — Ih escape width, and reinforces 
the above statement about the possible complexity of the GDR structure. The connection 
of these two different calculations supplies a value for the Tgdr that is close to the result 
obtained in the PDM calculations of the Ref. lll |. 

In short, our approach has been describing consistently the Ip—lh microscopic nature of 
the PDR excitation mode. The transitions in this region are composed by Ip— Ih excitations 
involving neutrons of 'skin'. The broad peaks are due to the coupling to continuum, which 
is extremely important in these fragile systems. Furthermore, we have exhibited an estimate 
about the direct escape of neutron from PDR. Regarding the GDR, due to small escape 
width obtained in our calculations, the results indicate that the more complicated excitations 
than Ip — Ih should also be important for the description of its microscopic structure. 
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